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Abstract 

We study the pairing correlations in a finite Fermi system from quantum entanglement point of 
view. We investigate the relation between the order parameter, which has been introduced recently 
to describe both finite and infinite superconductors, and the concurrence. For a proper definition of 
the concurrence, we argue that a possible generalization of spin flip transformation is time reversal 
operation. While for a system with indefinite number of particles concurrence is a good measure 
of entanglement, for a finite system it does not distinguish between normal and superconducting 
states. We propose that the expectation value of the radial operator for the total pseudospin can 
be used to identify entanglement of pairing. 
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Entanglement is a fundamental quantum mechanical propertyi which plays a central role 
in the quantum information theoryl. On the other hand, proposed measures of entanglement, 
including the entanglement of formation which quantifies the resources needed to create a 
given entangled statei, are generally not very proper for analytical calculations. Making use 
of a spin flip transformation, Wootter^l introduced so-called the concurrence to simplify the 
notion of entanglement of formation and Martin-Delgadoi extended and applied his results 
to a many-body problem, namely the BCS ground state of superconductivityi0. In this 
work, we examine the same problem for a finite system where the number of fermions is 
fixed. For this purpose, we make use of an order parameter proposed recently to describe 
both microcanonical and grandcanonical superconductors!. 

Experimental works on superconducting metallic islands at nanometer scale raised ques- 
tions about pairing correlationd^'B'S. For a bulk system, superconductivity can be described 
by a complex order parameter A. The equations have the symmetry that if A is a solution, 
then e*^A is also a solutionlli. However, in a finite system with fixed number of electrons, 
the order parameter A = (c_k|Ck|) vanishes since the operator does not conserve the num- 
ber of fermions. Here, c_ki and Ck^ are the annihilation operators for time reversed states 
I — k I) and | k f), respectively. In this case, superconductivity can be identified by non- 
vanishing number parity effect parameter Ap since the ground state energy of the system 
increases or decreases, depending upon whether the total number becomes odd or even, by 
addition of a new electron0'0. Recently, an order parameter has been proposed to unify 
the order parameter A of the bulk limit and the number parity effect parameter Ap of the 
nanoscopic superconductor^!. Using the pseudospin represent at ion0'El and the SU(2) phase 
stateJl^ a quantum phase has been defined for a superconductor with discrete energy levels 
along with modulus of the order parameter which becomes equal to Ap. As we go from 
the nanoscopic limit to the bulk superconductor it has been shown that the number parity 
effect parameter and the SU(2) phase go to the amplitude and the phase of the bulk order 
parameter, respectively. On the other hand, we can think of the long range order in the 
superconducting phase as an entangled state of Cooper pairs. In this paper we first discuss 
how to calculate concurrence as a measure of entanglement and we examine the relation 
between entanglement and the order parameter. We show that while for a system with 
indefinite number of fermions concurrence is a good measure of entanglement, for a finite 
system it does not identify the pairing correlations. As a possible solution we propose the 
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amplitude of our order parameter, which is nothing but the expectation value of the radial 
operator of the total pseudospin, to detect entanglement of Cooper pairs. 

For a finite Fermi system, such as a nanoscopic superconductor, energy levels are also 
finite and discrete and hence we can use a reduced form of the BCS mo del0 which was 
applied in nuclear physics and which has an exact solutionlll. The model Hamiltonian is 

^ = E ^cl-Cj- - ^ E CiTSV/iC/T (1) 

where g is the pairing coupling constant for the time-reversed states | j |) and \ j |), both 
having the energy ej. Here, c^^ (cja) is the creation (annihilation) operator for state | ja) 
where j G {!,... and a E {Tj i}- For the model Hamiltonian introduced above it has 
been shown that there exists a number parity effect, namely the ground state energy for 
even number of electrons is lower in comparison to neighboring odd number statesEi'00il 
including degenerate caseii. 

The key point in Wootters' formulation of concurrence is the spin flip transformation. 
For a pure state of a single qubit \ip) it is defined by 

1^) = ^y\r) (2) 

where {ip*) is obtained from lip) by taking complex conjugates of expansion coefficients and 
ay is Pauli spin matrix. For a single spin-| particle, this is nothing but the time reversal 
operation. The spin degree of freedom that we discuss here shouldn't be mixed with the 
pseudospin to be introduced below. For the many body case, a natural extension of the spin 
flip operation is the time reversal operation!. The action of the time reversal operator Ut 
on the creation operator is 

UtcIuI = (3) 
A similar relation holds for the annihilation operator. To find the transformed state 
(in the active picture) we can simply rewrite the transformed Hamiltonian (in the passive 
picture) and evaluate the corresponding eigenstate. 

To define the order parameter, we introduce the pseudospin variablesElS 

which obey the fundamental commutation relations of the SU(2) algebra 

±5,,sf. (5) 



Si , Sj 



It is possible to rewrite the model Hamiltonian as 

H = T.'^e,(s- + l)-9T.4sJ. (6) 

The mapping from the Fermi operators to the pseudospin operators is possible as long as 
all single particle states are doubly occupied. Since the original Hamiltonian (|l|) contains 
no terms which couple a singly occupied level to others, the only role of such states will be 
blocking from pairing interaction. Therefore, the summations in Eqn. are over doubly 
occupied and empty states only. Both the above mapping and the BCS wave functionil 
lack proper antisymmetrization due to separate treatment of singly occupied states, but 
since the model Hamiltonian (|l]) does not involve any scatterings into or out of such states, 
antisymmetrization with respect to interlevel pair exchange and intrapair electron exchange 
is sufficient. 

Given SU(2) algebra, for example the one generated by the components of the total 
pse.do.pin ope.ato. s ^ E. s„ we can int.odncji tKe .ad,a. operate, defined 

Sr = V S+S~ (7) 

and the exponential of the phase operator given by 

m=s 

E= \ S;sm + l){S;sm \ . (8) 

m=—s 

Here, | S; sm) is simultaneous eigenstate of and Sz operators with eigenvalues s{s + 1) 
and m, respectively. The label S has been introduced to distinguish them from the phase 
states to be defined below. For integer s or on the so called Bose sector, the eigenstate of E 
with eigenvalue exp(— z27m/(2s + 1)) is evaluated to be 
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0\ sn) = . exp 



27171 

i m 

2s + 1 . 



S; sm) (9) 



V2i+ 

and a similar expression holds for half integer s or in the Fermi sector. 

In terms of the radial and the exponential of the phase operators for the total pseudospin, 
it is possible to rewrite the interaction part of the Hamiltonian @ as —gs^EE'^Sr- Since E 
is unitary, we have EE'^ = I but we are going to keep E and E'^ without cancellation to 
introduce the phase properly. Now, we define {sr) and {E) as the amplitude and phase of the 
order parameter, respectively. It has been proven that (sr) becomes identical to the modulus 
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of the BCS order parameter in the bulk hmit while in the nanoscopic limit it reduces to the 
number parity effect parameter Apl. Furthermore, in the bulk limit, {E) becomes identical 
to the phase of the BCS order parameter. Next, we examine how the amplitude and the 
phase of the order parameter is transformed under the time reversal operation T. Equation 
(0) implies that the components of the total pseudospin operator s transform according to 

Uts^U^ = -s, (10) 

UxSyUrp = Sy 

UtsJJI = s^. 

The transformation has immediate consequences on the order parameter. First, the am- 
plitude (sr) remains unchanged. Second, since {E) = (s+)/(s^ ~ ^1 ^ Sz), the exponential 
of the phase expectation value acquires a minus sign or phase angle change by vr. This is 
consistent with the prediction of BCS mean field treatment!. Here, we note that pseudospin 
operator does not transform exactly like spin operator a whose components satisfy 

Ura.Ul = -a, (11) 

UxC^yUj^ = (Ty 
UtCTzUI = -Oz 



in the standard representation of Pauli matrices^. 

For a given state \^)^ the central quantity concurrence is defined byi 

C[^) =|(V^ I ^) I . (12) 

Since [s,-, s^] = 0, Sr leads to a good quantum number even for a finite system. The eigen- 
states and including the ground state of the model Hamiltonian will be of the form 

I ^ra) = ^c^\S\ sm) (13) 

s 

because the interaction term commutes with and while the single particle part com- 
mutes with the latter, only. In general, the total spin in is multiply degenerate. We can 
calculate the expectation value of the radial operator as 



(sr) = E I I' ^s{s + l)-m{m-l). (14) 



Since the problem is exactly solvable, Cg coefficients can be found numericallyE^. In terms 
of these coefficients we can write down the transformed state which will be same as (|13|) 
except that all coefficients will be replaced by their complex conjugates. Then we evaluate 
the concurrence as 

C^W=lEc'l- (15) 



The BCS ground state, which is superposition of states of the form (|T3D with different m 
values, corresponds to phase states | 9; sn) in our notation!. In other words, it is an extended 
state in m— space and C{BCS) can be calculated explicitly!. Here, we can evaluate the same 
quantity for the phase states. However, since the phase states are defined for a given s value, 
we need to generalize ( |T3| ) by 

\^e)^T.(^'"''\^m) (16) 

m 

from which we find that 

W.)«IEc'IIEe^""'l- (17) 

s m 

The second term implies that concurrence vanishes in contrast to the Fermi sea state. Hence, 
concurrence is a distinguishing parameter for entanglement of Cooper pairs. 

For a state with real expansion coefficients, assuming (|13D is normalized, concurrence is 
unity, i.e. it is same as unpaired state. This is the case for a system with fixed number 
of fermions. A simple and analytically solvable example is a system composed of a single, 
c?— fold degenerate energy leveliil. Therefore, microcanocial entanglement of pairing, which 
we define by following Martin-Delgado as the difference between concurrence values of the 
Fermi sea and the BCS ground state, vanishes. Although concurrence does not distinguish 
between the normal and the superconducting states, amplitude of the order parameter (s^) 
still identifies pairing correlations and hence it can be used as a signature of entanglement 
of Cooper pairs. 

In conclusion, for a superconducting system with indefinite number of particles, concur- 
rence vanishes while it is unity for the Fermi sea and hence it can be utilized to detect the 
existence of entanglement. On the other hand, for a finite system it is unity in the super- 
conducting state, too. However, the order parameter of superconductivity which we propose 
can still be used to identify entanglement of Cooper pairs. 
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